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Abstract 

We develop a general power counting scheme for the infrared limit of Landau gauge SU{N) Yang- 
Mills theory in arbitrary dimensions. Employing a skeleton expansion, we find that the infrared 
behavior is qualitatively independent of the spacetime dimension d. In the cases d = 2,3 and 4 
even the quantitative results for the infrared exponents of the vertices differ only slightly. Therefore, 
corresponding lattice simulations provide interesting qualitative information for the physical case. 
We furthermore find that the loop integrals depend only weakly on the numerical values of the IR 
exponents. 

Introduction: The infrared (IR) regime of Yang-Mills theory determines the fundamental non-pertur- 
bative properties of the non-Abehan gauge dynamics. In particular, the IR limit of Greens functions 
of the colored fields are likely connected to the problem of confinement as demonstrated within the 
scenarios of Kugo-Ojima [T] and Gribov-Zwanziger [21I3- Due to the intricate nature of the problem an 
investigation via all available methods is desirable. Such combined efforts via Dyson-Schwinger equations 
(DSE) [11 m m [7] , renormalization group (RG) techniques [8] and lattice gauge theory studies led 
during the last years to a coherent picture of the infrared regime in Landau gauge. Hereby, in the 
continuum approaches, the IR scaling Hmit of general Greens functions in four spacetime dimensions has 
been determined [lO]. It exhibits, as argued for in pj], a dominance of the gauge fixing part of the action, 
i.e. the Faddeev-Popov ghost. A simultaneous analysis of DSE and RG methods allowed furthermore to 
show that this IR fixpoint is unique [1.2]. 

The Greens functions of Yang-Mills theory in Landau gauge are studied extensively on the lattice, 
see e.g. p]. These analyses show the IR increase of the ghost dressing function required in both the 
Kugo-Ojima and the Gribov-Zwanziger scenario, as well as the IR scaHng of the gluon dressing function 
in accordance with DSE and RG analyses. The influence of finite volume corrections has been studied 
within DSE analyses on compact manifolds [13] which agree qualitatively with the finite IR limit of 
current lattice results. However, a confirmation that even the gluon propagator vanishes weakly in the IR, 
which is necessary in the Gribov-Zwanziger scenario, has not yet been possible in four-dimensional lattice 
simulations due to the large lattice sizes required to probe the IR regime. The situation is similar for the 
important but even harder determination of the scaHng behavior of vertex functions or the quantitative 
extraction of IR exponent. Lower dimensional lattices could allow to test qualitative aspects of these 
Greens functions in a substantially simpler setting. Naturally, such a surrogate study requires that the 
generic features are identical to those in four-dimensional spacetime. Recent lattice studies in two [H] 
and three spacetime dimensions [15] indeed suggest a power-law behavior similar to the four dimensional 
case. 

SU{N) Yang-Mills theory in three spacetime dimensions is also relevant from two other points of 
view. In the high temperature limit only the lowest Matsubara mode contributes, leading to a dimensional 
reduction and thus to an effective field theory in three dimensions. The qualitative aspects of this effective 
theory, in particular with regard to its confining nature, may be covered by ordinary three-dimensional 
gauge theory [16]. Similarly, Yang-Mills theory in the canonical quantization approach in Coulomb gauge 



[17] with an appropriate choice for the vacuum wave functional is structurally similar to the corresponding 
three-dimensional theory in Landau gauge. 

These results motivate to study the dependence of the infrared limit on the parameters of the theory in 
more detail. It will turn out that the infrared solution of Yang-Mills theory is determined by the spacetime 
dimension d and one parameter k that determines the IR exponents of propagators and vertices. In this 
work we establish a manifest power counting scheme for general vertex functions of Landau gauge Yang- 
Mills theory in arbitrary dimensions. We find qualitatively similar results as in the four-dimensional 
analysis given in [10] . Furthermore we discuss the dependence of the Dyson-Schwinger equations on the 
infrared scaling parameter k. In four dimensions it takes a value k « 0.595 within an approximation 
based on the propagator DSEs and using a bare ghost-gluon vertex. A possible dressing of this vertex 
could change this value slightly [6]. Therefore we study to what extent these dynamical building blocks 
in the DSEs depend on k, employing an analysis of the IR-dominant integrals. 



IR exponents for arbitrary d: In the following we will perform a scaling analysis for the IR regime 
of Yang-Mills theory in arbitrary dimensions. In contrast to the generic IR limit we study the limit 
where the coupling g is kept fixed and has no inherent scaling dependence. Instead in dimensions other 
than four it fixes the fundamental scale of the theory. This limit is more directly accessible in lattice 
simulations and has been analyzed in [HI [15] . A similar IR analysis has previously been performed for 
the physically important case (i = 4 in refs. [I0l[l2]. Although the scaling in arbitrary dimension may be 
abstracted from these results via the consideration of the appropriate canonical scaling dimensons, we 
chose to give a concise and nevertheless self-contained derivation of the general scaling relations and refer 
the reader to the corresponding work for details. We note, that our analysis is immediately applicable 
to the gauge sector of QCD since closed quark loops are IR-suppressed due to the finite current quark 
masses [18] and therefore there are no quark contributions to purely bosonic Greens functions in the IR 
regime. This has been verified in ref. [19] . 

Analog to the corresponding analysis in four dimensions [10] , the starting point for the IR analysis is the 
non-renormalization of the ghost-gluon vertex in Landau gauge [20] , which implies a finite vertex in the 
infrared [10] . This property depends purely on the transversality of the gluon propagator and is, as long 
as the ghost-gluon scattering kernel is not strongly IR-divergent, therefore vaHd in arbitrary dimensions. 
We come back to this point below. 




Figure 1: The DSE for the ghost propagator. 

The infrared behavior of the propagators and vertices well below its inherent scale (jr2/(4-<i) (respectively 
Aqcd in c? = 4) is determined via renormalization group arguments by scaling relations. The propagators 
of the gluons and ghosts 

are given in terms of dressing functions whose IR behavior is described by a power law ansatz 

Z(p2) = co,2 • , G{p') = C2,o • {p^Y'-° , (2) 

and similar for the vertices. Here we denote the IR exponent of a vertex with 2n ghost and m gluon 
legs by S2n,m and the corresponding coefficient by C2n,m- Whereas this coefficient is a constant for the 
propagators it is generally a function of 2n — m — 1 momentum ratios. 

For the integral on the right hand side of the ghost propagator Dyson-Schwinger equation, cf. fig. (H one 
can use the standard expression [6l [21] 



{27rY'^ ' ^ V / r(-i.i)r(-z.2)r(rf + j.i + z.2) 
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Figure 2: 1-loop part of the skeleton expansion of the DSE for the three-gluon vertex. 

which shows that it scales proportional to the external momentum. The left hand side of the ghost DSE, 
which consists only of the inverse dressed ghost propagator, is proportional to [p^)'^^-"^^ . The 1 comes 
from the canonical dimension of the ghost propagator. The scaling dimensions on the right hand side are 
d/2 from the integral, (5o,2 — 1 from the gluon propagator, — 1 from the ghost propagator, 1/2 from 
the bare ghost-gluon vertex, and 1/2 from the dressed ghost-gluon vertex which features no anomalous 
scaling. This yields as condition for the IR exponents 

1 - 62,0 = ^ + (5o.2 - 1 + 62,0 -^ + \ + \- (4) 

Defining the parameter k as k := —52,0, we recover in four dimensions the well-known result 5o,2 = 2k. 
In d dimensions we have O [6] 

<52,o = , (5o,2 = 2k + 2 - ^ . (5) 

A comparison between different dimensions via n is not possible directly, because n has different values 
for different d (SHU [16]. We discuss this point further below. 

In order to transform the infinite hierarchy of Dyson-Schwinger equations into a closed system we 
perform a skeleton expansion. This yields an infinite tower of graphs involving only primitively divergent 
vertices which we will analyze in a first step. The first order of the skeleton expansion of the DSE for the 
three-gluon vertex is depicted in fig. [21 We start with the ghost triangle which turns out to be one of the 
infrared leading diagrams. We have d/2 from the integral, 3(— k— 1) from the three ghost propagators and 
3/2 from the three ghost-gluon vertices. We subtract the canonical dimension 1/2 to get the anomalous 
IR exponent of the ghost triangle of the three-gluon vertex: 

<3'^ = ^ + 3(-«-l)+3i-i = -3^+^-2. (6) 

For the four-gluon vertex we can apply the same procedure. In fig. [3| we show the first order of its 
skeleton expansion. Here, also the ghost rectangle gives one of the IR dominant contributions. Simple 
counting of the powers yields for the IR exponent of this diagram 

5t^='^+A{-n-l)+A\ = -An+'^-2. (7) 

Now we are in the position to derive the IR exponents of the other diagrams of the three-gluon vertex 
DSE in fig. [21 The gluon triangle contains three propagators as well as two dressed and one bare three- 
gluon vertex, so including the loop integral and subtracting the canonical dimension the IR exponent is 
d/2 + 3(2k +1- d/2) -f 2(-3k -I- d/2 - 3/2) -I- 1/2 - 1/2 = 0. Correspondingly, the other corrections, 
including the two-loop diagrams that are not shown, are likewise subleading compared to the dominant 
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Figure 3: 1-loop part of the skeleton expansion of the DSE for the four-gluon vertex. 



ghost loop. The same holds for the four-gluon vertex. In general, the IR exponents of these purely gluonic 
diagrams are independent of the dimension and therefore identical to eqs. (13) and (15) in [lOj. 

With the scaling information of all primitively divergent Greens functions at hand we have all parts 
to calculate the IR exponent of an arbitrary n-point function. The different objects that can appear in 
corresponding general graphs are given in table [H In terms of these building blocks the IR exponent of 
an arbitrary vertex v is given by 

Sy =(/ - mj + Wo, 3 + "0,4)^ + (2mi - ~ 3wo,3 - 4t)o,4)K+ 

+ ^(2mj - 2ni + ^2,1 + 3 - 3uo,3 - 4wo,4 - 2u) , (8) 

where v is the canonical scaling dimension of the vertex. This formula can be simplified using standard 
relations between the number of loops, the number of vertices and the number of propagators. The 
number of loops is given in terms of the internal lines and vertices via the topological relation 

Z = + ni + 1 - {vis + i'o,3 + ^0,3 + t'S,4 + "^0,4) • (9) 

Further topological relations that connect the number of Hnes of a given species with the corresponding 
vertices allow to rewrite the expression in terms of the number of external lines which are denoted by m 
and n for gluons and ghost-antighost-pairs, respectively: 

m + 2mi = V2,\ + 3(^;o,3 + "0,3) + 4(wo^4 + wo,4) , (10) 
n + n, ^ V2,i ■ (11) 



object 


number of objects 


scaling dimension 


loop 


1 


d/2 


internal ghost line 


m 


S2.O - 1 = — K - 1 


internal gluon line 


mi 


(5o 2 - 1 = 2k + 1 - d/2 


ghost-gluon vertex 


V2.1 


1/2 


bare 3-gluon vertex 


^0,3 


1/2 


dressed 3-gluon vertex 


^^0,3 


So,3 + 1/2 = -3k + d/2 - 3/2 


bare 4-gluon vertex 


"0,4 





dressed 4-gluon vertex 


VOA 


60,4 = -4k + d/2 -2 



Table 1: Infrared behavior of the building blocks of general loop diagrams within the skeleton expansion. 
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The canonical dimension of the vertex is finally given in terms of its external lines 



4 — (2n + m) m 
v= L^2-n--. (12) 



Inserting eqs. JO]), |(l0]), ^ and ^ into eq. JH]) we get 

52n,rn = + 1 - ^0,4 ^ ^'0,3)^ + (4^0,4 + ^'"1,3 ^ "T- + n)n+ 

+ (2t>(;_4 + 2i;g 3 + 2n - 2) . (13) 

This formula is still dependent on bare vertices. However, when considering different terms of the skeleton 
expansion of an n-point function, one easily sees that terms containing bare three- and four-gluon vertices 
are not the IR-dominant ones because contributions from these terms are 

<3(3k-^ + 2), (14) 

<4(4'^-^ + 2), (15) 

and thereby always add a positive power of momenta. Consequently, we can neglect these diagrams, 
when counting only IR leading ones, i.e. we can set f § 3 = Uq 4 = in (fT3)) and obtain as a final result 
for the IR exponent of an n-point function the simple expression 

^2n,m = (ri - m)K + (1 - n) - 2^ . (16) 

This relation verifies a posteriori the assumption discussed above fig. [ij the ghost-gluon scattering kernel 
is exactly as IR-divergent as to make the ghost-gluon vertex finite in the IR, and this is independent of 
the value of the dimension d. 

It remains to discuss higher order corrections in the skeleton expansion. From eq. i|16p we can see 
immediately that there are infinitely many higher orders terms which have the same IR behavior as the 
leading term: Insertions which lead to higher orders only contribute with dressed vertices and propagators, 
and these were all considered in the derivation of eq. (flGl) . These insertions with their corresponding IR 
exponents are given by: 



T^^^^^^ (2K+l-|) + 2(-K-l) + 2i + f =0 

I I 

A 
B 



2(2«;+l-|) + (-K-l) + (-3«; + f -|) + i + | = 



W^^^^^ 3(2K+l_|) + 2(-3K-f f -f) + f = 



D 



2{2k + 1 - I) + (-4k + |-2) + f = 



5 



2,5 



2,0 



1,5 
1,0 
0,5 
0,0 
-0,5 



""0 1 2 3 4 5 6 

d 

Figure 4: The (mostly two) possible solutions for the IR exponent k as a function of dimension d, cf. 
also refs. [llj and [16]. The Hnear fits are exact in some ranges of dimension d and approximate the real 
solutions to good precision in others. The first branch contains the solutions k = 0.2, n = 0.3976... and 
K = 0.5953... for d = 2, 3 and 4, respectively, whereas the second one includes the solutions k = Q, k — 0.5 
and K = 1 [5]. 

It should be noted that in contrast to a uniform scaHng limit we have kept the (dimensionful) coupling g 
fixed. The IR power laws, eq. I|16p . are then valid for momenta (^4/(4-d)_ rjj^^ resulting momentum 

dependent couplings from the ghost-gluon vertex, a^^ , the three-gluon vertex, a^^ , and the four gluon 
vertex, a^" , (see [T^ for a definition) are then given by 

a<>\p') =f [Z,Ap')r G^ip') Zip") ^ (p2)(4-rf)/2 ^ (17) 
47r p^— ►o 

«4.(p2) (p2)] ^2(^2) (p2)(4-.)/2 (^g) 

47r p^^o 

Here Z2^i denotes the dressing function of the tree-level structure of the ghost-gluon vertex and Zo,3 
and Zo,4 the corresponding dressing functions of the three- and four-gluon vertices, respectively. To 
obtain these results the scaHng laws, eq. ifTB)) . have been used. From this one may better understand the 
usefulness of the skeleton expansion: apart from B all diagrammatic pieces given above contain dressing 
factors of the same multiplicity as the coupHngs, cf. eqs. l|17H19p . Thus all anomalous dimensions 
cancel and the canonical dimension (4 — d)/2 of these couplings is cancelled by the canonical dimensions 
and the extra loop. As a result the skeleton expansion indeed works independent of the value of the 
spacetime dimension d as could have been expected from the previously studied four-dimensional case 

Finally, let us compare the IR behavior of the Yang-Mills Greens functions in different dimensions. 
The corresponding values for the IR exponent k (using a bare ghost-gluon vertex) are displayed in fig. [H 
Apart from the region c? < 1 we always have two possible values, which arise from an IR analysis of the 
ghost and gluon propagator DSEs. We show the IR exponents of the ghost and gluon propagators and 
the three- and four-gluon vertices in table [21 Whereas the situation in four dimensions is not conclusive, 
yet, the solutions belonging to the first branch in fig.|4]have been found in lattice studies in two and three 
dimensions. These results for the IR exponents of the propagators, the ghost-gluon vertex ^22j and the 
three-gluon vertex in two dimensions |14j agree within errors with the corresponding values given in table 
[21 One can see that for this branch the qualitative behavior does not change in different dimensions. The 
other branch in fig. [H that starts at k = — 1 for d = and ends at k = 2 for d = 6 has not been seen in 
lattice simulations. We will give some additional arguments below that this branch may be unphysical. 



1 1 1 1 


1 ' 1 ' 1 ' 1 


" ■ Solutions for K 
- Fit: y = -1 + l/2*x 
— Fit: y = 215 - 2/5 *x 






— 


< 1 < 1 , 1 , 1 < 1 < 1 
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Dependence of the loop integrals on the IR-exponent: The Dyson-Schwinger equations involve 
loop integrals over bare and dressed vertex functions. As just discussed the corresponding dressing 
functions exhibit a power law with appropriate infrared exponents. The corresponding coefficients depend 
on the IR exponent k. For the respective coefficients of the propagators, eq. |[2|), which can be computed 
from the 2-point integral, eq. this K-dependence is shown in fig. [U where K-independent prefactors 
have been dropped. The physical requirement that the dressing functions have to be positive restricts the 
possible values of k to k < 0.5 in two, 0.25 < k < 0.75 in three and k > 0.5 in four dimensions ^23j. We 
note that all known solutions fulfill these requirements. The decisive point is, however, that although the 
curves in fig. [5] differ considerable, the physical values obtained from the Dyson-Schwinger solution are 
far away from these poles and thereby the K-dependence in their vicinity is rather mild and qualitatively 
similar in each case. 




Figure 5: The K-dependence of the unrenormalized dressing integrals appearing in the DSEs for the ghost 
(left) and gluon propagator (right). The dotted, dashed and solid lines show the curves for c? = 2,3 and 
4 respectively, see also ref. |[6]. The full and open dots represent the solutions of the first and second 
branch in fig. [H 



In general the vertex integrals feature more complicated tensor structures with an increasing number of 
independent tensor components. However, there are general methods to decompose such tensor integrals 
to standard scalar integrals [25 • The tensor integrals in the DSE for the ghost-gluon and three-gluon 
vertex reduce to 3-point integrals of the form 

where any additional momenta due to vertex functions in the numerator are included in the denominators 
for appropriate values Vi. A general expression for such scalar 1-loop 3-point integrals in arbitrary 
dimensions and with arbitrary powers of the propagators has been obtained in [21} I25j . An explicit 
evaluation of this result requires an expression for the resulting Appell function F4 whose defining series 
does not converge in the considered Euclidean regime. A converging expression can be obtained by 
analytic continuation and has been given in ^26j. Unfortunately, the result given there is sHghtly incorrect. 
The correct and lengthy expression which agrees with a direct numerical integration of these integrals 



Dimension 




4 


3 


2 


4 


3 


2 


K 




0.5953... « 0.6 


0.3976... « 0.4 


0.2 


1 


0.5 





Ghost 


-K-1 


-1.6 


-1.4 


-1.2 


-2 


-1.5 


-1 


Gluon 


2k+ 1 -d/2 


0.2 


0.3 


0.4 


1 


0.5 





3- gluon 


-3K + d/2 - 3/2 


-1.3 


-1.2 


-1.1 


-1.5 


-1.5 


-0.5 


4- gluon 


-AK + d/2-2 


-2.4 


-2.1 


-1.8 


-4 


-2.5 


-1 



Table 2: The dimension dependence of the infrared behavior of Yang-Mills Green functions. For each d 
the two solutions for k corresponding to the two alternative branches in fig. |4]are displayed. 
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will be given in a forthcoming publication [27j where we will also give the detailed results for the different 
tensor structures and their kinematic dependence. Here we will merely discuss the dependence of the 
overlap of the tree-level tensor with the IR-dominant ghost triangle. According to the discussed power 
counting, the ghost loop presents the IR-leading contribution to the Dyson Schwinger equation for the 
three-gluon vertex in fig. [21 Since this loop involves only ghost-gluon vertices that remain bare to leading 
order, its IR behavior can be analyzed semi-perturbatively using the scaling form of the dressed ghost 
propagators, see eq. We consider the special kinematic configuration given by the symmetric point 
— — r^. At this point the dressing function of the leading ghost loop correction to the three- 
gluon integral scales according to eq. ^ as Zq.s = co,3(k) • (p2-j-3K+rf/2-2 .^^j^gj-g ^j^g dependence of the 
coefficient on the scaling parameter k has been made explicit. The coefficient is shown as a function of k 
in fig. m where the full and open points represent the values for the first and second branch in fig. |4] and 
K-independent factors have been dropped again. 



0.1 



.05 



C0,3 



-0 .05 



-0 . 1 




0.2 



. 4 



. 6 



. 8 



Figure 6: The K-dependence of the overlap of the unrenormaHzed ghost triangle with the tree- level tensor 
at the symmetric point. The dotted, dashed and soHd lines show the curves for c? = 2, 3 and 4 respectively. 
The full and open dots represent the solutions of the first and second branch in fig. [H 



The vertex integral has zeros at k = in two, k = 0.5 in three and k = I'm four dimensions which coincide 
precisely with the solutions of the second branch. These zeros depend on the kinematic configuration and 
appear at higher values of k away from the symmetric point. (The corresponding scalar integrals, cf. eq. 
(|20l) . feature even poles at these values of k which precisely cancel for the tree- level tensor.) According 
to eq. (fTSl such zeros in the vertex dressing function would lead to an IR-vanishing coupling. Since 
Yang-Mills theory is apparently a strongly interacting theory, this indicates that these solutions might 
not be physically relevant. As in the case of the propagators the n dependence in the vicinity of the 
solutions of the first branch is similar in all dimensions. The negative value of the three-gluon vertex in 
two dimensions nevertheless leads to a positive IR Hmit of the coupHng a^^ in eq. lfT8|) . However, it would 
cause problems in its renormalization group fiow and might be an artifact of the skeleton expansion and 
cured by higher orders. 



Conclusions: We have studied generic features of the IR limit of SU (N) Yang-Mills theory established 
in [U [51 III [To] in more detail. As in the four-dimensional case, the IR behavior of Greens functions can be 
extracted via a skeleton expansion in arbitrary dimensions. We find that the IR Hmit of Greens functions is 
surprisingly insensitive on the spacetime dimension. This already constitutes the first important result of 
this study. As a consequence Yang-Mills theory in lower dimensions has a qualitatively similar IR limit as 
the four dimensional theory. Thus the confinement mechanism might possess identical features in different 
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dimensions. Corresponding lattice simulations should provide interesting qualitative information for the 
physical case. Indeed a recent study on large two dimensional lattices [H] finds infrared exponents in 
agreement with the value k = 0.2 , thus confirming the Gribov-Zwanziger scenario also in two dimensions. 

The results on the mild K-dependence of the remaining DSE solutions suggest that ghost dominance 
is a rather robust mechanism and should not depend on the details of the employed truncation scheme. 
This is further substantiated by the fact that the main non-linearities in the DSE system, which enable 
the non-trivial fixpoint, arise in the propagator equations whereas there is no non-Hnear feedback to 
leading order in the vertex equations. 
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